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Abstract 

Employing a perturbative QCD based parton cascade we calculate the elliptic flow V2 and its 
transverse momentum dependence V2(pt) for the gluon matter created in Au+Au collisions at 
v/ $nn = 200 GeV. To make comparisons with the experimental data at the BNL Relativistic 
Heavy Ion Collider (RHIC), parton-hadron duality is assumed. We find that whereas the integrated 
V2 matches the experimental data, the gluon (or pion) V2(pt) is about 20—50% smaller than the 
experimental data. Hadronization via gluon fragmentation and quark recombination seems to be 
the key to explaining the necessary jump of V2(pt) from the partonic to the hadronic phase. We 
also show that the elliptic flow values moderately depend on the chosen freezeout condition, which 
will thus constrain the shear viscosity to the entropy density ratio of the quark gluon plasma 
created at RHIC. 
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I. INTRODUCTION 



Elliptic flow is the key observable characterizing the collectivity |3j of the quark 

gluon plasma (QGP) that is potentially produced in ultrarelativistic heavy-ion collisions 
at the BNL Relativistic Heavy Ion Collider (RHIC) or/and at the CERN Large Hadron 
Collider (LHC). A fast transformation from the initial spatial anisotropy to the momentum 
anisotropy will indicate sufficiently strong interactions among quarks and gluons, which also 
provides the reason for a fast thermalization to a QGP with a small shear viscosity. The 
good agreements of the elliptic flow t> 2 data at RHIC J4, 5] with the results from calculations 
employing ideal hydrodynamics [6j suggest that the shear viscosity of the QGP created 
at RHIC is small Q. To find its lower and upper bound is an important but difficult 
task, because the elliptic flow is built up not only during the evolution of the QGP but 
also initially during the thermalization and finally, though possibly marginal, during the 
hadronic cascade before particles decouple kinetically. Moreover, the lack of understanding 
of the hadronization leads to an uncertainty in converting the partonic elliptic flow to the 
hadronic one that is measured. All these contributions to the final hadronic elliptic flow can, 
in principle, be investigated in a transport calculation including a partonic and a hadronic 
cascade and a hadronization mechanism 8j. Such studies will help to determine the true 



based parton cascade Boltzmann 



111 ], we have calculated the elliptic 



shear viscosity in the QGP phase. 

Recently, employing the perturbative QCD (pQCD _ 
Approach of MultiParton Scatterings (BAMPS) |9|, L 
flow t> 2 of the gluon matter produced in Au+Au collisions at the RHIC energy [12|. To com- 
pare our results with the experimental data the gluonic v 2 was converted to the hadronic one 
using the parton-hadron duality picture at a freezeout energy density of e c = 1 GeV fm~ 3 . 



Hence V2 lon (pr) 



gluon 



(Pt) at transverse momenta pt- Even though the hadronization was 



described in such a simple manner, it was found that when employing a QCD coupling con 



stant of a s = 0.6, the calculated t>2 matches the experimental data [13j, 1 141 ] . With a s = 0.3 
the results are about 20% smaller than the data. Hence, pQCD interactions can explain 

n n 

the large t> 2 buildup at RHIC. In addition, in Refs. [12| and [15] the shear viscosity to the 
entropy density ratio 77/s from the transport calculations was extracted and it was found 
that rj/s is between 0.15 for a s = 0.3 and 0.08 for a s = 0.6. These findings are similar to 
those obtained from viscous hydrodynamical calculations [la]. 
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This article provides details on our recent results on the elliptic flow from parton cascade 
calculations. In Sec. [Ill setups for numerical calculations are given. Section II I II shows the 
results for the transverse momentum dependence of the elliptic flow, ^(pt)) the transverse 
momentum spectra, the mean transverse momentum, and the final transverse energy per 
rapidity at midrapidity. The results from the BAMPS calculations are compared with the 
experimental data at RHIC assuming parton-hadron duality. Further possible improvements 
as well as a discussion about possible hadronization scenarios via gluon fragmentation and 
quark recombination are given. In Sec. IIVI the dependence of the elliptic flow on the cho- 
sen freezeout energy density is presented, which indicates the uncertainty of the T]/s ratio 
extracted from the BAMPS calculations. Finally we summarize in Sec. |V] In Appendix [AJ 
the calculation of the number of participating nucleons and the determination of centrality 
classes are given. 



II. BAMPS SETUPS 

BAMPS solves the Boltzmann equation for on-shell gluons with pQCD interactions, which 
include elastic scatterings and bremsstrahlung and its backreaction. The cross section of 
pQCD elastic scatterings is given by USD 



da" 9 9na 2 

" (qi + mlf ' () 

where denotes the perpendicular component of the momentum transfer in the center-o 



mass frame of the elastic collision. The interactions are screened by a Debye mass 0, 0, llSj , 



2 ( + \ a at f d3 P /( X '^P) w 
m D (x, t) = 7rd G a s N c J — — , (2) 

where d G = 16 is the gluon degeneracy factor for N c = 3. mp is calculated locally using 
the gluon density function /(x, t, p) obtained from the BAMPS simulation. The integrated 
cross section is 

2 Jo q± rfqi 2m 2 D (l + 4m 2 D /s)' K ) 

and s is the invariant mass of the collision. The factor 1/2 before the integral indicates that 
the outgoing particles in the collision are identical. 

The effective matrix element of pQCD inspired bremsstrahlung gg <-> ggg is taken in a 
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Gunion-Bertsch form [171 . 1 181 . [1 



\M 



9g 4 



aa^aaa I 



. - q ± ) 2 + m 2 D \ 



(4) 



2 (qi + m2,)»ki[(k 

where g 2 = 4ira s . h± and y denote the perpendicular component of the radiated gluon mo- 
mentum and its rapidity in the center-of-mass frame of the collision, respectively. k± = |kj_|, 
and A g is the gluon mean free path, which is calculated self-consistently The suppres- 
sion of the bremsstrahlung due to the Landau- Pomeranchuk-Migdal effect is effectively taken 
into account within the Bethe-Heitler regime using the step function in Eq. (j4]). Gluon ra- 
diations and absorptions are only allowed if the formation time of the process, typically 
t = coshy /k_L, is shorter than the mean free path of the radiated or absorbed gluon. Both 
the Debye screening mass and the gluon mean free path act as an infrared regulator. 
The total cross section for bremsstrahlung is obtained by the integral 

1 1 



gg^aaa 



-- J dr' 1 dr' 2 dr' 3 \M gg ^ ggg \ 2 (2n)^\ Pl + P2 - Pl - P ' 2 - P ' 3 ) 



(5) 



where dT\ = d 3 Pi /(27r) 3 /(2Ef i ). P \ and P 2 denote the four-momenta of the incoming gluons, 



and Pi 



1, 2, 3 denote the four-momenta of the outgoing gluons. For the backreactions 



we define a "cross section" by 



aaa^aa 



dr[dr' 2 \M 



aaa->aa I 



\2n)W 4 \ Pl + p 2 + n - p[ - p 2 ) 



(6) 



with \M ggg ^ gg \ 2 = \M gg ^ ggg \ 2 /d G . 

Interactions are simulated by the stochastic method [9|, |2(J, 
interaction (i.e., transition) probabilities 

°aa^aa ^ 
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22j, [23J according to the 



P 



aa-^aa 



aa^aaa 



aaa^aa 



VreV 



Vre\ 



Attest A 3 S 

a aa^aaa ^ 



test 



A 3 x 



l aaa-*aa 



At 



(7) 
(8) 
(9) 



8E 1 E 2 E 3 N 2 est (A3x) 2 ' 

for gg -> gg, gg -> ggg, and ggg -> gg processes, respectively. v rel = s/{2E l E 2 ) is the 
relative velocity between incoming gluons. A 3 x denotes the volume of the cell element and 
At the time step. The probabilities are calculated for every doublet and triplet in each spatial 
cell and are compared with random numbers between and 1 to determine whether or not 
interactions occur. If the random number is smaller than the probability for an interaction, 



the interaction occurs and new momenta of outgoing gluons will be sampled according to 
Eqs. ([1]) or (j4]). Details of the samplings can be found in Ref. If the random number is 
larger than the probability for an interaction, no interaction will occur. Particles that do not 
participate in interactions within At will propagate as a free streaming. These operations 
run over all the cells and time steps until the freezeout condition is reached. 

To reduce numerical artifacts, the cell length is set to be smaller than the gluon mean free 
path. For instance, in noncentral Au+Au collisions with an impact parameter of b = 8.6 fm, 
the transverse cell length is set to be a constant small value of Ax = Ay = 0.125 fm. The 
setup of the longitudinal cell length is refreshed before each new time step so that there are 
almost the same number of particles in each Az bin with Az = tftanh^^) - tanh(r/ s i)], where 
t is the time in the center-of-mass frame of a Au+Au collision and rj s = | ln[(i + z)/(t — z)\ 
denotes the space time rapidity. It turns out that the longitudinal cells are almost equidistant 
in r] s lOj, i.e., 1 77^2 ~ Vsi \ ~ const, which implies a nearly Bjorken-type space time evolution 
of the gluon matter. 

The accurate solution of the Boltzmann equation using the stochastic method according 
to Eqs. (J7j) — (J9]) is guaranteed if there are enough particles in each cell, because this condi- 
tion is assumed when deriving the interaction probabilities Eqs. (|7j) — ([9]) js]. Because the 
cell volume has to be small, the real gluon number in each cell is too small to fulfill the 
condition for using the stochastic method. We adopt the test particle method to amplify 
the particle number with a factor of N test . To leave the physical scale of the gluon mean free 



path unchanged, the interaction probabilities are accordingly reduced by N test and 



test ' 



respectively [see Eqs. (jZJ) — (|SJ)]. One "event" means N test events with parallel propagations. 
Particles from different events can interact but with the reduced probabilities, so that each 
particle still behaves as a true physical one. On the other hand, when calculating the real 
particle number and energy density, the values must be divided by N test . We set N test = 2400 
for b = 8.6 fm, which leads to \r] S 2 — f]si\ ~ 0.1 for the longitudinal cell length. 

In the present pQCD simulations the interactions of the gluons are stopped when the local 
energy density drops below e c . The value for e c is assumed to be the critical value for the 
occurrence of hadronization, below which parton dynamics is not valid. Because a realistic 
hadronization and a hadronic cascade are not yet included in BAMPS, a gluon, which ceases 
to interact, propagates freely and is regarded as a free pion according to the parton-hadron 
duality picture (5-function fragmentation of the gluon). Thus, e c determines the freezeout 
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condition. The implementation of a Cooper-Frye prescription for hadronization 24| and the 



subsequent UrQMD 



251 ] hadronic cascade after the QGP evolution is in progress. 



Because the implementation of Bose enhancement into transport calculations is still tech- 
nically difficult, gluons are treated as Boltzmann particles for simplicity. If the gluon system 
is in thermal equilibrium, the energy density is e = 48T 4 /7r 2 . In Ref. 12J we have chosen 
e c = 1 GeV fm -3 , which leads to a critical temperature of T c = 200 MeV. To see the pos- 
sible dependence of the elliptic flow on the freezeout density we choose, in this article, an 
additional moderately lower value of e c = 0.6 GeV fm -3 that gives T c = 175 MeV. 



The initial gluon distribution 
Ref. 12J: an ensemble of minijets 



br BAMPS calculations is the same as that chosen in 



26j with transverse momenta greater than p = 1.4 GeV, 



produced via semihard nucleon-nucleon collisions. The value of po is chosen by matching the 
final transverse energy per rapidity from the parton cascade calculations to the experimental 
data [121 ] . We use Glauber geometry with a Woods-Saxon profile and assume independent 
binary nucleon-nucleon collisions. The inclusion of quarks in BAMPS is straightforward and 
is in progress. We do not expect large changes in the integrated elliptic flow t>2, because 
the quark amount at RHIC is initially about 20% only. As quarks interact weaker than 
gluons, we expect that Vi will be slightly smaller than that obtained in the present pure 
gluon matter approach. More discussion about this issue appears at the end of the next 
section. 



III. ELLIPTIC FLOW AND DISCUSSIONS ON HADRONIZATION 
A. Results 

The BAMPS calculations for Au+Au collisions at the RHIC energy y/s^N = 200 GeV 
are carried out for a set of discrete impact parameters from b = 2 to 11 fm (see Appendix 
|A|) . The corresponding number of participating nucleons A part (6) is calculated within the 



wounded nucleon model 27|. The details are presented in Appendix |A] The results from the 
BAMPS calculations with the given impact parameters are used to obtain the elliptic flow 
as a function of the transverse momentum ^(pt) and the transverse momentum spectra in 
different centrality classes. Calculations with randomly sampled impact parameters are not 
appropriate for BAMPS, because a single event is already extremely time consuming for 
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large test particle numbers N test . 

Figure [I] shows the elliptic flow t> 2 = ((pi — v Vllv V) as a function of the number of 
participating nucleons, iV part . The points are STAR [13| and PHOBOS [14j data for charged 
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FIG. 1: (Color online) Elliptic flow vs iVp art for Au+Au collisions at ^s^n = 200 GeV. The 
points are STAR [J and PHOBOS [1J| data for charged hadrons within \rj\ < 0.5 and \r]\ < 1, 
respectively, whereas the curves with symbols are results for gluons within |?7| < 1, obtained from 
the BAMPS calculations with a s = 0.3 and 0.6 and with two freezeout energy densities, e c = 0.6 
and 1 GeV fm" 3 . 



hadrons within the pseudorapidity intervals \r]\ < 0.5 and |?7| < 1, respectively. The symbols, 
which are connected with colored straight lines, are results for gluons within \r]\ < 1 (rj being 
identical to the momentum rapidity y for massless gluons) from the BAMPS calculations for 
two values of the coupling constant a s and for two values of the critical energy density e c 
for the freezeout. Especially, the blue curve with open circles and the red curve with open 
squares are the results for (a s = 0.3, e c = 1 GeV fm -3 ) and (a s = 0.6, e c = 1 GeV fm~ 3 ), 
respectively. These are exactly the same as shown in Fig. 2 of Ref. [Pj. The new results 
for e c = 0.6 GeV fm -3 will be discussed in the next section. 

Figure [2] shows the elliptic flow as a function of the transverse momentum ^(pt) for the 
most central 50% of Au+Au collisions at ^snn = 200 GeV. The points are PHOBOS data 



14j for charged hadrons within < r\ < 1.5, whereas the curves with symbols are results 



for gluons within \r)\ < 1.5 from BAMPS calculations. Figure [3] shows ^(pt) m different 
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FIG. 2: (Color online) Elliptic flow vs transverse momentum for the most central 50% of Au+Au 
collisions at ^/snn = 200 GeV. The points are PHOBOS data [14] for charged hadrons within 
< rj < 1.5, whereas the curves with symbols are results for gluons within < 1.5, obtained from 
the BAMPS calculations with a s = 0.3 and 0.6 and with e c = 0.6 and 1 GeV fm -3 . 

centrality classes for particles within \r)\ < 0.5. The PHENIX j^, 28] and STAR data [l^] 
are depicted for comparisons with the results from the BAMPS calculations. The latter are 
plotted until pr = 4 GeV, because the statistical errors at higher pr are too large. 

Although the integrated v 2 of gluons from the BAMPS calculations with a s = 0.6 and 
e c = 1 GeV fm -3 agree perfectly with the experimental data (see the red curve with open 
squares in Fig. [T|, the differential v 2 (pt) (the red curve with open squares in Fig. [2] and red 
curves in Fig. [3]) are 20—50% smaller than the data, especially at intermediate momenta 
1.5 GeV < px < 4 GeV. Correspondingly, the px spectra of gluons, shown in Fig. HI are 
flatter than those of hadrons at RHIC. The larger the centrality, the larger is the difference 
between the data and the results from the calculations. 

Discrepancy with the RHIC data is also seen in the mean transverse momentum, (pr), 
which is shown in the upper panel of Fig. 0as a function of A^ part . The points are STAR data 
{30} within < 77 < 1, whereas the curves with symbols are the BAMPS results for gluons 
within \r]\ < 1. The mean p? of gluons is larger than that of charged hadrons. On the other 
hand, the total transverse energy of gluons per rapidity at midrapidity, shown in the lower 
panel of Fig. [5j agrees with the experimental data including hadronic and electromagnetic 
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FIG. 3: (Color online) Elliptic flow vs transverse momentum for different centrality bins of Au+Au 



2^1 and STAR data for charged 



collisions at ^/sWn = 200 GeV. The points are PHENIX [4|, 
hadrons within < 0.35 and \r]\ < 0.5, whereas the curves are results for gluons within \r/\ < 0.5, 
obtained from the BAMPS calculations with a s = 0.3 and 0.6 and with e c = 0.6 and 1 GeV fm~ 3 . 

components [30]. We notice that the results of dEx/dr] with a larger coupling constant of 
a s = 0.6 are 7—15% smaller than those with a s = 0.3 (also smaller than the data) because 
of the larger longitudinal work during a less viscous hydrodynamical expansion. However, if 
quarks are included in BAMPS calculations, the transverse energy will be slightly enlarged. 



B. Discussion 



dEx/dri agrees rather well with the data but (pt) is larger than the data. This thus 
indicates that the gluon number per rapidity is smaller than that of the hadron number 
per rapidity. The parton-hadron duality with one gluon to one pion seems too simple to 
describe the complex hadronization. To obtain (px) that agrees with the data, each gluon 
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FIG. 4: (Color online) Transverse momentum spectra for different centrality bins of Au+Au col- 
lisions at y / J/v r /v = 200 GeV. The points are PHOBOS data [23] for charged hadrons within the 
rapidity interval 0.2 < y n ± < 1.4, whereas the curves are results for gluons within |y| < 1.5, ob- 
tained from the BAMPS calculations with a s = 0.6 and e c = 1 GeV fm -3 . To compare with the 
experimental data, the gluon number is reduced by a factor of 3. 

is expected to "fragment" to 1.5—2 pions on average. The larger the centrality (peripheral 
collisions), the more fragments are needed, because initial hard gluon jets are less quenched 
in peripheral collisions because of their small reaction size. 

Possible fragmentation of each gluon to more pions suggests ^^(pr) ~ vf non (npT) and 
dNP ion (p T ) ~ dm lnon (np T ) with n = 1.5-2. The obtained vf on (p T ) for low p T < 1 GeV 
would come closer to the experimental data. Because the integrated V2 is determined mostly 
by the lower pj- particles, the integrated v-i of pions will be nearly the same as that of 
gluons. On the other hand, the V2 lon (pT) for p? > 1 GeV would be still 20—50% smaller 
than the data. The discrepancy is larger for more peripheral collisions (see Fig. EJ). It 



seems that quark recombination models 3l|, [32j, |33|, |34j for hadronization are appropriate 



for intermediate pt gluons, if one assumes that each gluon must first be converted to a 
quark- ant iquark pair g —>■ qq via t>2 Uark (j>r /2) ~ wf luon (pr). Then the quark recombination 
models would give i>2 lon (pr) ~ 2t>2 Uark (p T /2)) ps 2ff luon (p T ). 

We also must stress that the results discussed above are obtained from the present calcu- 
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FIG. 5: (Color online) Mean transverse momentum vs iVp ar t (upper panel) and transverse energy 
per pseudorapidity vs -/Vp art (lower panel) in Au+Au collisions at ^fsmq = 200 GeV. The points 



are STAR data [30j] within < rj < 1, whereas the curves with symbols are results for gluons within 
1^1 < 1, obtained from the BAMPS calculations with a s = 0.3 and 0.6 and with e c = 0.6 and 1 
GeV fm" 3 . 

lations without quark degrees of freedom. The inclusion of quarks will certainly change the 
freezeout temperature, because the freezeout energy density is proportional to the degener- 
acy factor of partons, e c ~ (da + dq)T^, where do = 16 for gluons and g?q = 24 for quarks 
with two flavors. For the same freezeout energy density the freezeout temperature will be 
a factor of 1.26 smaller when the quarks are included and if all the partons are chemically 
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equilibrated. Although the initial energy is slightly enlarged when including quarks, we do 
not expect much stronger transverse flow at freezeout, as long as the shear viscosity of the 
parton system is the same with or without quarks. Therefore, the inclusion of quarks might 
reduce the mean px by a factor of 1.26, which then would lead to a fair agreement between 
the calculated and the experimentally measured mean px via parton-hadron duality. If the 
integrated vi is still comparable with the experimental data, then V2{pt) will also come to 
a closer agreement with the data. Work in this direction is in progress. 

The parton-hadron duality picture is consistent with the hadronization in a thermalized 
expanding matter, which assumes that both constituents of the bulk matter at the beginning 
and at the end of the hadronization will be distributed according to thermal statistics. If 
the temperature does not change much during the phase transition, the integrated v 2 , v 2 (px) 
and the pt spectra of pions are almost the same as those of partons. Although the energy 
density of pions is smaller than that of partons, the matter will expand continuously during 
the hadronization, which leaves the total energy, and the total particle number and, thus, the 
entropy of the constituents approximately unchanged. Within the picture of a thermalized 
expanding source the dependence of ^(pr) on the hadron mass was also understood [(J and 
agrees with the RHIC data at lower px < 1-5 GeV 



On the other hand, the STAR data 
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32 



33 



35] confirmed quark recombination models [3l|, 



34J ] for hadronization of quarks with intermediate px ~ 1—2 GeV. These models lead 



to V2 S,dron (pr) ~ nv2 Uark (pr I 'n) for n = 2 for mesons and for n = 3 for baryons, which is 
obviously different from the picture of a thermalized expanding source. How the transition 
from one to another hadronization scenario proceeds is not clear. Also, how gluons convert 
into quarks and antiquarks in the recombination models has to be understood and modeled. 

We emphasize that how the partons hadronize is essential to obtain a complete picture 
of the generation of elliptic flow and to understand the transformation of the v 2 {px)^ the pt 
spectra, and the mean px from the partonic to the hadronic phase in a consistent manner. 
Quantitative calculations modeling hadronization may give a more detailed understan ding 
of this subject. As a starting point, an algorithm employing a Cooper-Frye prescription [24j 
will be implemented in BAMPS in the near future. Detailed comparisons with results from 
recent calculations using viscous hydrodynamics 36| will be given. 

Finally, the gluons in the present BAMPS calculations are assumed as Boltzmann par- 
ticles rather than bosons because of numerical difficulties. Bose enhancement increases the 
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gluon density by a factor of £(3) ~ 1.202 for a thermal gluon system. Moreover, the en- 
hancement at low pt is considerably larger than that at large pt- This may also have a 
moderate effect on the dependence of the elliptic flow and thus will be investigated in 
more detail. 



IV. DEPENDENCE OF THE ELLIPTIC FLOW ON THE FREEZEOUT CONDI- 
TION 



The elliptic flow i> 2 will be continuously built up if the spatial anisotropy of the particle 
system still exists and sufficiently strong interactions among particles are present. The 
buildup of t>2 does not necessarily mean that v% becomes always larger. After the time when 
the eccentricity crosses zero the system will turn spatially anisotropic again but with the 
sign opposite to the eccentricity of the initial one. Further evolution will then decrease v^. 
Whether or not a generation of negative occurs depends on whether the kinetic freezeout 
of the particle system happens after or before the time when the initial eccentricity crosses 
zero. 

In this section we discuss the dependence of the elliptic flow on the freezeout condition in 
Au+Au collisions at the RHIC energy. In the present BAMPS calculations the kinetic and 
chemical freezeouts are assumed to occur at the same time when the local energy density 
drops below the chosen cutoff value of e c . The cutoff e c serves as a critical energy density, at 
which gluons are converted to hadrons. After the phase transition no hadronic interactions 
take place in the present version of BAMPS. We have already speculated that a realistic 
modeling of hadronization may change the pattern of ^(Pt) i n various ways. Moreover, a 



hadronic cascade simulating the decoupling stage should be included |37|, |38|, |39( . With the 
hadronic cascade, the kinetic freezeout happens later than the hadronization and the final 
elliptic flow will likely be different from that obtained in the present BAMPS calculations. 
However, this difference is expected to be rather marginal, because the viscosity in the 
hadron gas is much larger than that in the QGP. 

Although the freezeout is considered in a simple way in the present BAMPS, it is useful 
to study the uncertainty in the calculated t>2, if the value of the critical energy density e c is 
varied to change the lifetime of the QGP. If the freezeout condition constrains the elliptic 
flow, then it also constrains the to be deduced shear viscosity of the QGP. For this purpose, 
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we have performed calculations with e c = 0.6 GeV fm -3 in addition to those with e c — I 
GeV fm -3 . The results are shown in Figs. [U-13] and [5] by green curves with open triangles 
(a s = 0.3) and purple curves with open diamonds (a s = 0.6). 

We realize that the v 2 results with a s = 0.3 and e c = 0.6 GeV fm -3 (green curves with 
open triangles) are almost identical to those with a s = 0.6 and e c = 1 GeV fm -3 (red curves 
with open squares). Stronger interactions or longer QGP phase leads to the same final values 
of v 2- Figure [6] shows the v 2 generation as a function of time in a Au+Au collision with 
an impact parameter of b = 8.6 fm. No decrease of v 2 is observed. This indicates that the 




b=8.6 fm 

a =0.3, e =1 GeV f m 3 

S C 

a =0.6, e c =1 GeV fm" 3 

a =0.3, e =0.6 GeV fm" 3 

S ' C 

- a =0.6, e =0.6 GeVfm" 3 



t (fm/c) 

FIG. 6: (Color online) Time evolution of the elliptic flow from the BAMPS calculations in Au+Au 
collisions at y / ijv r /v = 200 GeV at an impact parameter of b = 8.6 fm. 



freezeout occurs before the initial spatial anisotropy vanishes. The saturation of v 2 begins 
at 2.5 fm/c for e c = 1 GeV fm -3 and later at 3 fm/c for e c = 0.6 GeV fm -3 . The continuous 
increase of v 2 after 2.5 fm/c in the case for e c = 0.6 GeV fm -3 is as strong as that before. 
The difference of the pressure gradient between the x and y directions is still large at 2.5 
fm/c, whereas at this time the freezeout at e c = 1 GeV fm -3 is nearly complete. Therefore, 
the uncertainty in the final elliptic flow due to the different freezeout condition is not small. 
In addition to the hadronization, the time scale when the hadronization occurs affects also 
the absolute value of v 2 as well as the transverse momentum dependence of v 2 . 

A detailed study of the elliptic flow is important, because this collective effect of QCD 
matter quantitatively constrains the shear viscosity of the medium. In Ref. 12| we have 
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demonstrated within the BAMPS calculations that the shear viscosity to the entropy density 
ratio rj/s is nearly constant in time and approximately depends only on the coupling constant 
a s : T]/s ~ 0.15 for a s = 0.3 and rj/s « 0.08 for a s = 0.6. From Fig. [T]we see that with 
either a s = 0.3 and e c = 0.6 GeV fm~ or a s = 0.6 and e c = 1 GeV fin" 3 the same values of 
the final v 2 are obtained and which agree very well with the experimental data. Therefore, 
according to the present study, rj/s is most probably lying between 0.15 and 0.08. 

If the freezeout occurs at a lower energy density (or temperature), one would expect that 
the transverse momentum spectra will become softer and the final transverse energy will be 
smaller. From the lower panel of Fig. [5] we see that dEr/dr] for e c = 0.6 GeV fm -3 are only 
slightly smaller than those for e c = 1 GeV fm -3 . The differences in the transverse momentum 
spectra (not shown) are also tiny. The further decrease of the local energy density from 1 
to 0.6 GeV fm -3 due to the longitudinal work done by the pressure is marginal, because 
at the late stage of expansion the system becomes dilute and thus the work done is small. 
Free streaming increases the transverse flow and thus effectively decreases the local energy 
density. Therefore, a lower energy density cutoff for the freezeout does not lead to a much 
smaller total transverse energy and mean transverse momentum, as seen in Fig. [5j On the 
other hand, the ^(pr) for ce s = 0.6 and e c = 0.6 GeV fm -3 , shown by the purple curve with 
open diamonds in Fig. [2] and by purple curves in Fig. [3l agree well with the data at the low 
Pt < 1.5 GeV, whereas the integrated t>2 (see the purple curve with open diamonds in Fig. 
[[]) overestimate the data. 

We note that within the same description for the gluon interactions the energy loss of 
high pt gluons is found to be in good agreement with the results from the GLV formalism 
when a s = 0.3 is used [do}. Hence both jet-quenching phenomena and the buildup of elliptic 
flow are rather well described in BAMPS parton cascade calculations. 

V. SUMMARY 

Employing the pQCD based parton cascade BAMPS that includes pQCD bremsstrahlung 
and its backreaction we have calculated the elliptic flow v 2 in Au+Au collisions at ^/Jnn — 
200 GeV. Hadronization is assumed to happen at a critical energy density, e c . A gluon 
is then fragmented into a pion according to the parton-hadron duality assumption. The 
hadron w 2 is therefore identical to the final gluon v 2 within the present BAMPS calculations. 
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We found that whereas the final gluon t> 2 from the calculations with a s = 0.6 and e c = 
1 GeV fm -3 agree well with the data, the transverse momentum dependence i> 2 (pr) are 
20—50% lower then the data in each centrality class. The final gluon transverse momentum 
spectra are harder than the data and also the final gluon mean transverse momenta are 
40—100% larger than the data. However, the final gluon transverse energy per rapidity at 
midrapidity agrees with the data. This indicates that the hadronization process via the 
parton-hadron duality is not justified. A realistic hadronization that transforms a gluon to 
1.5—2 pions on average would be expected for low p T gluons, whereas at intermediate p?, 
quark recombination would be a better scenario of hadronization. In addition, the inclusion 
of quark degrees of freedom would reduce the final parton mean px and would lead to a better 
agreement between the calculated and measured i> 2 (pt)- This must be demonstrated in a 
new version of BAMPS calculations including the quark dynamics and employing various 
hadronization scenarios. 

The value of the final gluon t> 2 depends on the freezeout condition, i.e., e c . The pressure 
gradient difference in the x and y directions at e c = 1 GeV fm~ 3 is still large enough to 
further increase the elliptic flow if e c is changed to a smaller but still reasonable value. We 
observed that the final gluon t> 2 and w 2 (pt) are almost the same in the calculations with 
either a s = 0.6 and e c = 1 GeV fm~ 3 or a s = 0.3 and e c = 0.6 GeV fm~ 3 . Stronger 
interactions or a later freezeout leads to the same elliptic flow. This outlines the uncertainty 
in the extraction of the shear viscosity in the QGP: The shear viscosity to the entropy 
density ratio will most probably be between 0.08 and 0.15. 

We furthermore conclude that adding quark degrees of freedom into the dynamical evo- 
lution of the QCD matter with a detailed understanding of the hadronization of gluons and 
quarks will be helpful in explaining the viscous facets of the final hadron elliptic flow and in 
extracting the shear viscosity to the entropy density ratio of the QGP. The results presented 
in this article motivate more detailed investigations of this issue in future works. 
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APPENDIX A: NUMBER OF PARTICIPATING NUCLEONS iV part AND CEN- 
TRALITY CLASSES 

The number of participating nucleons N pait (b) in an A+B collision at an impact parameter 
of 6 is calculated within the description of wounded nucleons [27| : 



N paIt (b) = / d 2 sn paxt (s,b) , (Al) 

where 

n part (<f, 6) = T A {s) [l - e -^T B (M] + TB (b - g) [l - e -*HT A (3)] (A2) 

on denotes the nucleon-nucleon total inelastic cross section with diffraction production ex- 
cluded and is set to be oh = 42 mb. T A (s), also Tg(s), is the thickness function defined 

as 

/oo 
dzp A (s,z), (A3) 
-oo 

where pa(s,z) is the single nucleon density and / d 3 r pA^r) = J d 2 sT A (s) = A. We use the 
Woods-Saxon function for the nucleon density 

PA (f)=PA(r)= n l RA , (A4) 
1 + 

where d = 0.54 fm, n = 0.17 fm~ 3 , and R A = 1.12A 1 / 3 - 0.86A- 1 / 3 = 6.37 fm for A = 197 
of a Au nucleus. 

We use N part (b) to make a relation between the intervals of the centrality class and the 
impact parameter b. For an interval of b, [61562]) the average iV part and the average impact 
parameter can be calculated by 

ft dbbN pSLIt (b) 2 

dbb °2~ °i Jb i 

ltdbb 2 2 6 3 - 6 3 

{b)l[hM = ^^ = sW^M- (A6) 

The number of events, Advents; within [61; 62] is proportional to b 2 , — b\. Therefore, 

-^events I [6 i; 6 2 ] _ 6| - b\ 

N** ent8 Max{6j +1 - 6|} ' 1 ] 
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where N^ ents is the maximum of N events within all the intervals [bj] bj + i\, j = 1, 2, • • •. 

We obtain the intervals of b corresponding with each centrality class by tuning bjs, so that 
(N pai t)\[b d -b j+1 ], and Events I [6 j ;6 3 - +1 ]/^ent s from Ec l s - (EIJ) — (EJ) are comparable 

with the experimental data given for each centrality classes. Tables |T] and [III show the 
correspondences of the intervals of b with the centra 
the calculated values and the data from PHOBOS 



ity classes and the comparisons between 
29 1 and STAR Q. 



TABLE I: Correspondence of the impact parameter intervals with the centrality classes by match- 
ing the PHOBOS data Q. 



Centrality 



[bj]b j+1 ] (fm) 



0-6% 6-15% 15-25% 25-35% 35-45% 45-50% 
[0; 3.16] [3.16; 5.23] [5.23; 6.95] [6.95; 8.31] [8.31; 9.46] [9.46; 9.98] 



(N paxt ) PHOBOS [29 
(A part > Eq. (HU 



AWents/A^ n ts PHOBOS [29] 
N cvcnts /N^ cnts Eq. (H2l) 



344 ± 11 276 ±9 200 ± 8 138 ± 6 93 ± 5 65 ± 4 

350 282 208 147 101 74 

0.505 0.883 0.997 0.997 1.0 0.5 

0.477 0.829 1.0 0.99 0.98 0.48 



TABLE II: Correspondence of the impact parameter intervals with the centrality classes by match- 
ing the STAR data Q. 



Centrality 
[bj]b j+1 ] (fm) 


0-5% 
[0; 3.22] 


5-10% 
[3.22; 4.68] 


10-20% 
[4.68; 6.59] 


20-30% 
[6.59; 8] 


30-40% 
[8; 9.36] 


40-50% 50-60% 
[9.36; 10.42] [10.42; 11.07] 


{N paxt ) STAR [13] 


352 ±6 


298 ± 10 


232 ± 10 


165 ± 12 


114 ±12 


75 ±11 


46 ±9 


(A part ) Eq. (TX5D 


349 


293 


226 


160 


108 


69 


46 


(b) (fm) STAR [13] 


2.3 ±0.2 


4.2 ±0.3 


5.9 ±0.3 


7.6 ±0.4 9.0 ±0.5 


10.2 ±0.5 


11.3 ±0.6 


(b) (fm) Eq. flMD 


2.15 


3.99 


5.69 


7.32 


8.70 


9.90 


10.75 



BAMPS calculations are carried out for a set of discrete impact parameters 6 = 2, 2.8, 
3.4, 4, 4.5, 5, 5.6, 6.3, 7, 7.8, 8.6, 9.6, 10.4, and 11 fm. The average of an observable O in 
one centrality class is calculated by the integral 

f£dbb<D(b) 



<0)k 



Sgdbb 



(A8) 
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using the trapezoid formula. 
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